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In this paper we give a first principles microphysics derivation of the nonequihbrium forces between 
an atom, treated as a three (iiinensiorial harmonic oscillator, and a bulk dielectric medium modeled as 
a continuous lattice of oscillators coupled to a reservoir. We assume no direct interaction between the 
atom and the medium but there exist mutual influences transmitted via a common electromagnetic 
field. By employing concepts and techniques of open quantum systems we introduce coarse-graining 
to the physical variables - the medium, the quantum field and the atom's internal degrees of freedom, 
in that order to extract their averaged effects from the lowest tier progressively to the top tier. The 
first tier of coarse-graining provides the averaged effect of the medium upon the field, quantified by a 
complex permittivity (in the frequency domain) describing the response of the dielectric to the field 
in addition to its back action on the field through a stochastic forcing term. The last tier of coarse- 
graining over the atom's internal degrees of freedom results in an equation of motion for the atom's 
center of mass from which we can derive the force on the atom. Our nonequilibrium formulation 
provides a fully dynamical description of the atom's motion including back action effects from all 
other relevant variables concerned. In the long-time limit we recover the known results for the 
atom-dielectric force when the combined system is in equilibrium or in a nonequilibrium stationary 
state. 
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I. INTRODUCTION 

It is a remarkable fact that the dominant force between neutral atoms at short distances (e.g. the London force 
between two hydrogen atoms dominates at distances < 1mm) arises from quantum fluctuations. Indeed, at zero 
temperature for distances larger than the wavelength associated with an atom's first optical resonance the interaction 
originates from vacuum fluctuations. The last decade saw intensified research in this area of fluctuation-generated 
forces between atoms and between an atom and a conducting or dielectric surface. This was brought about by increased 
high-precision capability in the manipulation of trapped atoms in cavities and optical lattices [U , superconductivity 
experiments [3] and the design and operation of nano- and micro-electromechanical devices [IHB], amongst others. 
These advances which hold the promise of ushering in a new era of quantum engineering also made possible a wide 
range of theoretical inquires, such as the measurements of non-Newtonian forces [Q n-IL . utilization of experimental 
systems out of thermal equilibrium [12l HI] , and the investigation of quantum field theoretical effects in tabletop 
experiments. The expanded capability of ultrafast, high intensity lasers and high-precision manipulation techniques 
on atoms make it possible to observe real time processes which in turn demands a new theoretical framework to treat 
systems far from equilibrium. In our research program which began with the two earlier papers [161 El] and continues 
with the present work we apply the methods of nonequilibrium quantum field theory [T4| within the open quantum 
system |15| conceptual framework to tackle these problems. This method can provide a fully dynamical description of 
(non-stationary) systems far from equilibrium under the influence of various environments, or acted upon by different 
noises, going beyond the traditional mean field or linear response treatments. 

The state-of-the-art as far as we can discern from the literature is the so-called macroscopic quantum electrody- 
namics (MQED) which describes electromagnetic field fluctuations in a lossy medium. In order for the system to 
remain in thermal equilibrium the energy absorbed by the dielectric medium is compensated for by a stochastic 
forcing term which is added by hand in a manner consistent with the fluctuation-dissipation relation (FDR) |20) . 
Previous works employing MQED (or its variants) [ 2T| - [24l [27| [29H3T| (to name a few) have been skillfully employed 
to study atom-surface forces (or spontaneous emission) for a plethora of experimental setups. The key limitation of 
this technique is that it requires the system to be in a steady-state, or at least in local thermal equilibrium in order to 
apply the FDR. The method we use reproduces these earlier results for systems under equilibrium or nonequilibrium 
steady-state conditions, but is capable of treating fully nonequilibrium conditions including the effects of relaxation, 
and arbitrary atomic motion including full back-action from the field and the medium in a self-consistent manner. The 
FDR emerges from our model under steady-state conditions as a consequence of our energy conserving microscopic 
formulation in distinction to being required as in the formulation of Rytov's theory |25) . 

The techniques of MQED have been adapted to describe nonequilibrium atom-surface forces for two scenarios. In 
[26j a full treatment of the interactions between the body-assisted field in equilibrium and an atom in an arbitrary 
state was undertaken. There, it is the time-devolpment of the atomic state that characterizes the nonequilibrium 
behavior. In distinction Antezza et. al. |27j have successfully applied a nonequilibrium generalization of MQED 
to describe atom-surface forces when the field and surface are not in global thermodynamic equilibrium yet remain 
stationary (see [55] for a unified treatment of both scenarios) . This generalization relies crucially on the validity of a 
'local source hypothesis', which assumes no spatial correlation of the fluctuating polarization that drives the field in 
MQED. This ad hoc hypothesis is equivalent to ignoring interactions among the micro constituents of the dielectric 
medium. When the temperature of the body is much higher than the interaction energy among the medium's micro 
elements we believe that the 'local source hypothesis' should be an excellent approximation, and that the permittivity 
can safely be assumed to be local in space. But at low temperatures or when the coherence length of fluctuations in 
the medium become large this approximation is expected to break down and new techniques are needed to probe the 
fully nonequilbrium regimes. 

A key challenge in this endeavor is to understand the effects of dissipation on the field. Introducing a complex 
permittivity to the theory by hand provides dissipation but the energy of the field is not conserved, and so the field 
cannot be quantized in terms of energy eigenmodes. A first step was taken in [331 134j who considered quantization of 
the field in a dielectric half-space for the case of a real and frequency-independent permittivity. For real permittivity 
there is no dissipation and the field can readily be quantized but the permittivity violates the Kramers-Kronig relations 
[551 and so engenders acausal response. 

Our ultimate aim in this paper is a first-principles derivation of the atom-surface force based on the microphysics 
of the dielectric medium, the quantum field, and the atom. At the micro-level the total combined system is energy 
conserving. Dissipative effects arise at the macroscopic level after coarse-graining the detailed information in one 
or more subsystems. We employ the model of Huttner and Barnett [37] who described the micro-elements of a 
dielectric material by a continuous lattice of harmonic oscillators (from now on we'll refer to it as the matter field) 
coupled to a reservoir (see Appendix A of [38] for the matter described in terms of a fermionic theory). The use 
of a reservoir provides a microscopic method for introducing dissipation into the remaining degrees of freedom in 
the total system. For such a model composed of field + dielectric + reservoir the problem can be solved exactly by 
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Fano's diagonalization f39' when the coupHng between respective components is bihnear. In distinction to that work, 
we are not particularly interested in the microscopic details of the dielectric but only need to capture the averaged 
effect of the medium on the quantum field and the atom. By invoking the concepts and techniques for open quantum 
systems we coarse-grain the medium by tracing over the dielectric variables leading to a complex permittivity (in the 
frequency domain) which accounts for the dielectric's response to the field, and the absorption and emission of energy. 
In addition, the fluctuations of the oscillators which make up the dielectric manifest as a stochastic polarization that 
drives the field, and when the system is in equilibrium serve to balance the dissipative losses. From this microscopic 
viewpoint we see that a complex permittivity, quantifying material response to an external field, cannot be added 
to the theory freely unless the fluctuations of the same degrees of freedom are accounted for. It can be shown that 
at this level the semi-classical equations of motion for the field under the influence of the medium at late times ( 37 ) 
take the exact same form as in MQED (Lifshitz) theory when the medium is assumed to be in a thermal state. The 
microscopic approach we take offers a unique vantage point to see that the stochastic polarization driving the field put 
in by hand in MQED without field quantization actually arises from what is equivalent to the media's fluctuations. 
One could interpret the field fluctuations in MQED as being induced by a fictitious matter field. 

The approach we adopt is similar to that of 40; where a path integral formulation was used to derive the effective 
action describing the medium- influenced dynamics of the electromagnetic field. For the specific case of a dielectric 
half-space our results can be compared with [?T] who generalized the results of Carniglia and Mandel to frequency- 
dependent and lossy permittivities. They calculated exactly the dielectric -t- field dynamics using the Wiener-Hopf 
method and a sum over diagrams. We go beyond these results by providing a fully nonequilibrium treatment and 
apply these results to the atom-surface force. 

The paper is organized as follows. In Sec. 2 we describe the microscopic model. In Sec. 3 we find the medium 
altered equations of motion for the quantum field. From there the permittivity of the medium is identified as well as 
a stochastic force which accounts for the response and fluctuations of the dielectric. In Sec. 4 and Sec. 5 successive 
layers of coarse-graining are performed to obtain the dynamics of the atom's trajectory. In Sec. 6. we consider 
the specific case of the atom-surface force in a dielectric half-space and compare with the known results from the 
literature. We adopt the Einstein summation convention and natural units throughout fi. = c = fcs = 1. 



II. MICROSCOPIC MODEL 

The action describing the entire system S\z, Q, A^^, P, X^] is the sum of eight terms: 



S[z, Q, A^ P, X] ^Sz + Sq + Se + Sm + Sx + St'^ + + , (1) 

with five free actions pertaining to the five dynamical variables and three interaction actions. Here 1) Sz is the action 
for the motion of the atom's center of mass with coordinate z and total mass M under the influence of an external 
potential F[z(A)]: 



dt 



(2) 



where t is the atom's worldline parameter. 2) The internal degrees of freedom of the atom are modeled by a three 
dimensional harmonic oscillator with coordinate Q and natural frequency D, with action 



Sq[Q] = f 



dt 



Q'{t)~n'Q\t) 



(3) 



where ji is the oscillator's reduced mass. 3) The dynamics of the free photon field is described by Se where 'E' stands 
for the electromagnetic field. 



SEW]^\jd^xF^^''F^,. (4) 

where F^^ = d^j^A^ — d^A^ is the field strength tensor, with A^ the photon field's vector potential and / d^cc = 
//.^ dt J d^x. 4) The matter may be described by a continuous lattice of harmonic oscillators with natural frequency 
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uj which is meant to model the polarization of the medium, the coordinate of each oscillator is described by the vector 
field, P, with action 



Sm[P] = \j^ d^x[P\x) - u:^P\x)\. (5) 



The subscript V on the integration denotes that the spatial integration is restricted to the volume containing the 
matter. 

5) Each oscillator comprising the matter field is coupled to a reservoir. The reservoir is composed of a collection of 
oscillators at each point with frequency dependent mass I{u) and coordinates X^, with natural frequency u 

Sx[X.] = \j^d^xj dv I{v)[kl{x) - u^X\x)]. (6) 

The interaction of these parties is specified by the three interaction actions. 6) The interaction between the internal 
degree of freedom (dof ) of the atom and the field is 

Q,z\=qJdX Q'{X)Ei{X, z{X)) (7) 

where q represents the electronic charge. For the case of the matter the dipole moment of each oscillator is coupled 
with the local electric field with a coupling Qp 



5,^f[A^p] = ^p / d^x pUx)e,{x). 

Jv 



(8) 



For our model, which considers only the coupling of the electric field to the local polarization of the matter, there is 
no magnetic response for the medium and so the permeability can assume its vacuum value, /Xo, throughout. It should 
be noted that we have not included interactions among the elements of the dielectric which will provide a spatially 
local form for the permittivity in the macroscopic Maxwell's equations. 

Finally, each oscillator composing the matter is coupled to a reservoir with frequency dependent 'charge' g{v) which 
will provide dissipation and noise 



(9) 



SfJ[P,X.]= I d'x r dv g{v)Pi{x)Xi{x). 
Jv Jo 

III. FIELD EQUATIONS IN THE PRESENCE OF A MEDIUM 

In this section we will show how coarse-graining over the medium degrees of freedom leads to a permittivity and a 
classical stochastic source responsible for an additional induced component of field fiuctuations. Here, when we refer 
to medium we mean the combined reservoir + matter system. For now let us forget about the atom and field and 
focus entirely upon the matter and reservoir. 

A. Reservoir-Reduced Density Matrix 

Consider the time evolution of the density matrix describing the medium (matter + reservoir system). f){t) = 
U{t,ti)p{ti)U^t,ti) from some initial state at time <j to a final time t. By considering the matrix elements in an 
appropriate basis we can express p{t) as a product of path integrals 

p{Pf, P;, X^f, XI f, t) = P'''' VP exp{«(5M[P] - Sm[P'\ + P\ - s^n^W'.P'])} 

JCTP 

^ .x,.,,x„, ^^^{i^Sx[XA - Sx[Xl] + SfJ[P, X.] - SfJ[[P', Xim (10) 

., J OTP 
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where the density matrix depends upon the electromagnetic field through the interaction Sl^^ and CTP stands for 
closed-time path. We will assume that the initial state of the total system factorizes allowing the convenient notation 

r'^f'^f ^ r ^ r ^ r'^f ^ r^f ^ ^ ^ 

/ VY= VY, VYl VY VY' p{Y,,Yl:U) (11) 

JCTP J J JYi JYf 

adopted to keep long expressions compact where the integrals over Yi and Yl trace over all configurations of the field 
Y at the initial time. To capture the averaged effect of the reservoir on the matter we trace over it's final variables 
leading to the reservoir-reduced influence functional, Fx [P, P'] 

Trx{p{Pf,P'f, X,f, Xlf, t)} = r^''' VP exp{z(5M[P] - Sm[P'] + S[/,[A>^,P] - Cf K', ^'])} 

JCTP 

xFx[P.P']. (12) 

The reservoir-reduced influence action is given by 



i\nFx[P,P']^Sfp[P,P'] = 



d^x' 



P'- {t, x)e^''ix, x')P+{t',x) + -^P'-{t, x)&"{x, X')P- [f, x) 



(13) 



where <&'"'\x,x') and &"{x,x') are the retarded and Hadamard Green's functions for the reservoir, respectively, and 
the superscripts -I- and — denote semi-sum and difference variables deflned by P"*" = {P + P')/2 and P~ — P — P' 
respectively. Here the unprimed and primed quantities denote the forward and backward histories in the closed-time- 
path (or Schwinger-Keldysh) sense 

The reservoir kernels take the explicit form 



and 



(5''''\x,x')= diy {g^{iy)/iyI{iy))sm:y{t-t')9{t-t')S^{x- x') x &V 
Jo 



&"{x,x')^ I dv {g^{v)/vI{v))co\h{l3xvl2)cosv{t^t')5^{x- x') x (^V 
Jo 



(14) 



(15) 



where we have assumed that the initial state of the reservoir is thermal with inverse temperature /3x- The matter 
is not spatially correlated because the reservoir oscillators at one position do not interact with their neighbors. The 
retarded Green's function represents the response of the reservoir to an external held and so corresponds to its 
susceptibility. The Hadamard function, which can be derived from the symmetric two-point function of the reservoir's 
variables, quantifles reservoir fluctuations. This can be seen by noting that when the time arguments are coincident 
the Hadamard function becomes 



&"{t,t)^ / diy{g\,y)/,yI{iy)){AXUt)) 
Jo 



(16) 



where (&^{x,x') = &^ {t,t')5^{x — x') and X^{t) is the displacement of the r^th reservoir oscillator at a point. 

By the assumption that the reservoir is initially in thermal equilibrium these kernels satisfy the fluctuation- 
dissipation relation. Taking the Fourier transform of the time dependent component of each kernel we find 



e-''\w,x,x')= / di. (5'(^.)M(i.)) 
Jo 



+ i— [ 5{uj — v) — 5{uj + v) 



5^{x-x') 



(17) 



and 



©^(a;,f,f') = / dv {g^{v)/vI{v))coi\i(Pxv/2) 
Jo 



5{uj — v) + 6{u] + v) 



(53(f-x') 



(18) 
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noting that the imaginary part of the reservoir susceptibiUty ©''^*(a;, x, x'), which describes the absorption of energy 
by the reservoir, is balanced by the reservoir's fluctuations 

&"{oj,x,x') = 2coth(/3xw/2)Im[©''^*(a;,f,x')]. (19) 

To see explicitly how reservoir fluctuations influence the matter we appeal to a Gaussian path integral identity first 
suggested by Feynman and Vernon [44] . Note the complex modulus of the influence functional can be written as 

= J VUnix]exp!^i J^d^x ^x{x)-p-{x)^ (20) 
where 'P[^j(:] takes the form (with a normalization constant) 

V[ix]^e^V^~ j^d^x j^d^x' S,k e^(x)©^i(a^,a^')ex(^')} (21) 
allowing the reservoir-reduced influence functional to be expressed as 

^x[P,P'] = J V^xV[^x]exp!^iJ^d''x p-ix)[(,\,{x) + J^d^x'i&^'\x,x')P'+{x')]y (22) 

In this process we have replaced the kernel describing quantum and thermal fluctuations in the reservoir with the 
classical variable S,x which drives the field in the same manner as an external stochastic force. To retrieve information 
about the reservoir's fluctuations it is necessary to integrate over the functional distribution ^[^x] which is positive 
definite because the kernel <8jj^ is symmetric and positive. Therefore we interpret the field £,x as a classical stochastic 
force or noise, driving the matter with probability distribution described by 'P[(,x]- Due to the Gaussianity of 
all of its moments are specified by the mean ($.x) = and the variance ({^xi^)'^xi^')}) = 6^''<Sh{x,x') where 

{...)^^^JV^xP[ix]i...). 

Now that we have coarse-grained over the reservoir we can derive the semi-classical equation of motion for the 
matter field, ignoring for now its interaction with the field, to see the reservoir's averaged effect on the matter. A 



saddle point approximation of ( 12 ) after the fluctuation kernel has been replaced by a stochastic forcing term gives 



+uj'^P^ - j^' dt'(S''''\t,t')P^{t') ^Cxit)- (23) 

The effect of the reservoir is now manifest; first, the addition of the integral kernel in the equation of motion leads 
to dissipation in the free oscillations of the matter field, and second, the fluctuations of the reservoir drive the 
matter field through the source term on the right. For the specific case of Ohmic reservoir spectral density, that is 
g'^{v)/vl{v) — 2^v/'K, we find 

pk ~2pk ^ ^pfe(t') ^ (24) 

where the integral kernel leads to a dissipative term and a frequency renormalization cj^ = cj^ — 2jS{0) and the 
fluctuation kernel defining the distribution of (*) becomes 



<QH(t,t') = — [ du vco\h{l3xy/'2.)cosv{t~t')5^{x-x') feK 
7^ Jo 



B. Medium- Reduced Density Matrix 



(25) 



Consider now the time evolution of the density matrix elements describing the combined field -I- medium system, 
where the reservoir has already been integrated out, expressed as a product of path integrals 
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X / I^Cx^Kx] / VPcxp^i 

JCTP 



JCTP 



Sm[P] " Sm[P'] + SZ W\P] - Sr^, [A^ ,P']+ I d'x ■ P- 

V 



(26) 



The action 5m contains the dissipation term from the reservoir-reduced influence action 



1 



Sm[P]^Sm[P] + ;, I d^x / d^x' P\x)@-''\x,x')P,ix'). 
^ Jv Jv 



(27) 



Upon tracing the density matrix of the combined field + medium system over the matter variables we obtain a reduced 
density matrix p^. — Trpp that accounts for the averaged effect the medium has on the field. For linear coupling 
between the medium and the field and an initially Gaussian matter state the path integrals can be evaluated exactly 
yielding the medium-reduced density matrix |45j . The result of tracing over matter variables gives 



pM'},A'}-t) 



(28) 



CTP 



where the medium-reduced influence functional J-mIA'^ , A^^ ] accounts for the averaged effect of the medium on the 
field. The influence action Sfp, which relates to the influence functional as J^m = J T^£,x'P[£,x] '^^p{iSfp}, is given by 



Sfp[A>',A'']^ / d'^x / d^x' Er{x)\f'*{x,x') 



nlE'+{x')+npCxix') 



+ -nj,~g"{x,x')E^-ix'] 



(29) 



The kernels g'''^*(a;, x') and {x, x') are the matter's retarded and Hadamard functions, respectively. In this case the 
retarded Green's function does not follow directly from the commutator of the matter's polarization operators because 
the semi-classical equation of motion for the reservoir- influenced polarization field is not self-adjoint. It can be solved 
for classically or by considering the commutator of the polarization operator with its adjoint with respect to (23 1. 



Because the matter -|- reservoir system is linear the Heisenberg equation of motion for the polarization operators is 
the same as the semi-classical equation (24). The symmetric two-point function of the homogeneous solution of this 
equation subject to the appropriate boundary conditions provides g^{x,x'). Because the oscillators comprising the 
matter do not interact and the medium is assumed to be isotropic, the polarization fluctuations within the medium 
are not spatially correlated i.e. g{x,x') = g{t,t')S^{x — x'). Under these conditions the influence action simplifies 
greatly 



]^ d^x dt' E-{t,x)\ ~gret{t,t') 



nlE^+{t',x) + np^^x{t,x) 



+ -nl~gHit,t')E^-{t',x)l (30) 



'V Jti 

where the explicit form for gret{t,t') and gH{t,t') are given below for Ohmic spectral density of the reservoir 



9ret{t, t') = le-V2(t-t') _ ty^t _ i') 



(31) 



,-7/2(t+t'-2t,) 



3th(/3pa;/2) 



Co^ cosuj[t — t') 



— cosuj(t + t' - 2ti) + -ujs\\iu:{t + t' - 2U) 
4 2 



(32) 



where UJ = ^o)^ — 7^/4 and ti is the initial time at which the reservoir and matter begin to interact. Note in particular 
that the Hadamard function is not stationary and thus describes a nonequilibrium process. It can easily be verified 
that when the matter-reservoir coupling goes to zero these kernels give back the expected free forms. 

As with the case of coarse-graining over the reservoir described in the previous section we see that tracing over the 
matter degrees of freedom leads to two terms that are nonlocal in time: the term containing the retarded Green's 
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function characterizes the response of the dielectric to the field {gret plays the role of the susceptibility), and the 
one containing the Hadamard function quantifies the effects that quantum and thermal fluctuations (/3p being the 
matter's inverse temperature) of the polarization within the medium has upon the field. 

As with the case of the reservoir we can represent the fluctuation kernel in terms of a stochastic force. This allows 
us to write the complex modulus of the influence functional as 



where T'[Cm] takes the following form modulo a normalization constant 



V[S,m] = exp 



d'^x / d'^x' S,jCm{x)9h 
V Jv 



which defines (Uif(a;), CSfl^c')}) - S^^'-gnit^nS^x - x'). 
Using 7^[^Af] the medium-reduced influence functional becomes 



(33) 



(34) 



+ dt' r'M{^pE^+{t,x)+ex{t,x)) 



(35) 



Now we turn to the stochastic effective action S^[A^^ , A'^ ] for the field which is the sum of the free- field action and 
the influence action after the terms quantifying the fluctuations of the medium have been replaced with stochastic 
forces. As we work with the photon field in the path integral a gauge-fixing prescription must be adopted in order 
to prevent summing over gauge equivalent orbits. The Fadeev-Popov trick could be employed but the ghost fields 
introduced do not couple to the electromagnetic field in flat space and only contribute an overall multiplicative 
constant. Furthermore, the currents in our microscopic model which couple to the photon field are conserved. So we 
are free to choose any gauge we wish to evaluate the path integral [47j. As a result the Green's functions which appear 
will be gauge dependent. However, as a consequence of current conservation any gauge choice will give equally valid 
descriptions of physical processes. With this freedom we can express the stochastic effective action for the field, after 
choosing the temporal gauge (A° = 0). At this level we can make a comparison with MQED by allowing tf and —ti 
to go to infinity, assuming the field evolves to a steady-state and taking the Fourier transform to find 



1 

2^ 



duj' d^x 



- (V X A~{-uj\x)) ■ (V X A+{uj',x)) 



+A i-Lj', x) ■ A+{uj', x)uj'^{l + nlgretiu}', x)) + VtpE (-w', x) ■ gret{uj', x)^x{uj', x) 



(36) 



where in the infinite time limit the matter's stochastic force is exponentially suppressed 0. Note that we've 

added the explicit position dependence of the medium's retarded Green's function because its support is restricted to 
the volume containing the dielectric. 

The stochastic semi-classical equations of motion for the field can be derived from the saddle point condition for 



the medium-reduced density ( 28 ) matrix which gives 



V X V X A{uj', x) — uj' {1 + flpgreti^', x))A{uj', x) = iVlpijj'greti^' : a;)Cx(w', x). 



(37) 



In this form the role of the coarse-grained medium is evident. By comparing with the macroscopic Maxwell's equations 
we identify the permittivity as e{ui',x) = 1-1- fl^greti^' , x) and observe that the fluctuations of the medium drive 
the field through the stochastic force(s) ^x{x) {^m{x) drives the field in addition to the reservoir when the system is 
not in a steady-state). For Ohmic spectral density of the reservoir the frequency dependent permittivity corresponds 
with the Lorentz oscillator model 



e{uj' , x) 



uj'{uj' + i"f) — up- 



x^V. 



(38) 
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When the restoring force of the matter vanishes we obtain the Drude model and the plasma model when 7 — )■ 
in addition. In this form it's clear that the matter-field coupling flp can be interpreted as the medium's plasma 
frequency. 

One can see from equation (37) a striking similarity with MQED. Due to the linearity of our theory this is not 
surprising. Indeed if one were to proceed from this point treating the field semi-classically and choosing all space 
to be filled with a dielectric material, albeit in vacuum this dielectric is fictitious, one would exactly reproduce the 
predictions of MQED using (371 and choosing the dielectric to be in a thermal state. However, it is important to 
note that in this case the stochastic field ^-f^ represents the fluctuations of the medium only and does not include the 
intrinsic fluctuations of the field. After the next level of coarse-graining described in the following section the intrinsic 
quantum fluctuations of the field will enter which is different from the induced fluctuations from interaction with the 
dielectric medium. 



IV. FIELD-REDUCED DENSITY MATRIX 



In the last section we showed how the medium influences the field, and for the particular case of local fluctuations 
we found that the stochastic semi-classical action for the field takes the same form as the noninteracting field (with 
frequency dependent velocity) driven by an external current (stochastic force). 

A detailed description of the field's microstate is unimportant in the description of the atom-surface force. By 
coarse-graining over the medium-influenced field we can incorporate the averaged effect of the field on the atom's 
trajectory without a specific knowledge of the final field state, leading to the reduced density matrix for the atom 



Pr{zf,z'f,Qf,Qy,tf)^ Vz VQeM^iSz[A"Sz[A + SQ[Q]-SQ[Q']} 

JCTP JCTP 

/fA'^ A'^ 
VAI / ' ' ' VAf^ eMSfJ [z, Q,A'^]- S^,^ ,Q\A^']) }Tm K , ] (39) 
JcTP 

where the first integral in the second line traces over the final field configurations. 

For the assumed initially factorized state the path integrals over the field can be evaluated exactly if the initial state is 
Gaussian yielding the field-reduced influence functional, J^^ [ J^, J'" , ^] = J VS^x'Piix] J 'D^m'P[£.m] exp{iSfp [J'', ,£,]} 
expressed in terms of the influence action, Sfp , given by 



b^;^{x,x'){r+{x') - ntcix')) + -b'^Ax,x')r-{x') 



(40) 



where £,{x) = ^m{x) + J*;' dt' gret{t, t')^x{t', x). The current density J^^ in (40 1 comes from the atom-field interaction 
and takes the explicit form 



Jt^(x) ^^q dX Q'{X)k'^5\x°' - z"(A)) 



(41) 



where the derivative operator = — 9o'7f + diTj^ yields the electric fleld when contracted with the vector potential 
E'^ = K^^A^ and also enforces current conservation d^K,'^f{x) = {—dodi + dido)f{x) = 0. The integral kernels Z)^^* 

and -D^ are the retarded Green's and Hadamard function for the medium-altered electromagnetic field which result 
from solving the semi-classical equations of motion (for the retarded Green's function sourced by a delta function). 
The retarded Green's function describes the classical electrodynamical propagation of the field in the presence of the 
dielectric material, and the Hadamard function describes the field's intrinsic quantum fluctuations. 

In the temporal gauge the semi-classical equation of motion for the field's retarded Green's function in the presence 
of a dielectric medium satisfies 



e^\e"^\dad^bi;,l\x,x') + —bl^\x,x') + - 



dh gret{t,t2;x)—b:^\t2,x,x') = S,kS\x - x') (42) 



where Habc is the Levi-Civita symbol (Roman indices refer to s pati al components). The solution to (42) gives the 
particular solution, A^ , to the semiclassical equation of motion (37) for the electromagnetic field 
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A^{x) = - / d''x'b]l\x,x')^e{x'). (43) 

As was noted previously the Heisenberg equations of motion for the field operators take the same form as the 
semi-classical equation of motion because of the linear coupling 

e"\e'""&5aa„,i„(a;) + ^i,(a;) + — dt2gret{t,t2;x)—Mt2,x)^—^,{x) (44) 



therefore we can use the homogeneous solution to (44) to construct the symmetric two-point function of the field 
operators to give the Hadamard function 

&^ix,x')^({Alix),Aiix')}). (45) 
Finally, the field-reduced density matrix takes the form 



JcTP JcTP 

V. ATOM'S INTERNAL DOF (Q)-REDUCED DENSITY MATRIX 

At this point we have the density matrix that describes the dynamics of the atom's trajectory and its internal degrees 
of freedom under the infiuence of the medium-altered field. As with the previous parties it is only the averaged effect 
and not the microscopic details of the atom's internal degree of freedom which we need for the description of the force. 

As our last tier of coarse-graining we trace over Qf the oscillator's internal dof to obtain the oscillator-reduced 
density matrix which characterizes the dynamics of the atom's trajectory determined by its interaction with all 
remaining parties 

pz{zf,z'f;tf)= r^^'' Vze'^^-^'^-^-^^^^FzizJ]. (47) 

JCTP 

The environmental influences are now packaged in the oscillator-reduced influence functional !Fz z"] with their 
back-action accounted for in a self-consistent manner. The coarse-grainings are summarized in Table |T] at each tier 



TABLE I. Summary of Coarse-Graining and Detailed Physics Remaining at Each Tier 



Tier Coarse-Graining Influence Functional Detailed Physics Remaining Variables 

all parties z,Q,A^,P,X^ 

1 Trx J^x atom -I- field -)- matter z.Q^A'^.P 

2 Trp Tm atom -I- field z,Q,A^ 

3 Tta J'i atom z, Q 

4 Trg Tz atom's motion 2 



the influence of all lower tiers on the remaining degrees of freedom is packaged in the influence functional. 

We proceed by evaluating Tz\z^ z!\ perturbatively to lowest order in the atom-fleld coupling. In the fully dynamical 
case an exact calculation is not possible as backaction of the field on the dynamics of the atom's dipole moment enters 
as a third time derivative and includes multiple refiections between the dielectric medium and the atom. At leading 
order in powers of the atom-field coupling, an expression for the force can be found which neglects the radiation 
reaction to the atom's dipole moment and the effects of multiple reflections. 

To ease the notational burden we define 
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JcTP 



(48) 



which represents the noninteracting time-dependent expectation value with respect to the oscillator's initial state. 
With this simplification the oscillator-reduced influence functional can be compactly expressed in terms of quantum 
and stochastic expectation values 



^tSip[z,z'] de_f 



Tz[z, z] = 



iSf 



(49) 



which introduces the influence action, Sip. 



point is determined by the equation 



A saddle point approximation of (47 1 gives the semi-classical equation of motion for the trajectory where the saddle 



bz^ (r) 



ScGEA {z-, z' 



M-ik{T)+dkV[z{T)] = Ut), 



(50) 



the coarse-grained effective action, Scg ea [z, z!\ , is defined as Sz \z\ — Sz [^] + Sif [z, z!\ , and the influence force, (r) , 
is given by 



bz^-ir) 



S 



IF [z,z 



= fkir). 



(51) 



We now wish to evaluate the influence force perturbatively for an expansion in terms of small atom-field coupling. 
We begin by expanding both sides of (49 1 . Because Sfp contains terms linear and quadratic in the atom- field coupling 
we find to O(g^), 



1 



J^z[z, z'] = l + iSip + 0{q') « 1 + z {(Sfp)^) - - (((5i=.)2) J + 0{q') 



(52) 



where the subscript ^ in (...)^ represents that both the and stochastic expectation values have been taken. On 
the right hand side there are two contributions: the term linear in Sfp reduces to 



{{SFf)J^ = J d'^x dV^(j^-(x)J''+(x'))„ + - {J^~{x)r-{x'))^b^,{x,x') 

while the term linear in ^-^ vanishes. The quadratic term in Sfp reduces to 



(53) 



{{{SfFf)o), = / d^^' A dy {j^-{x).r-{y))^D'^;;{xy)Dl^^{v,y')K'in^ {Cix'Wiy')) 



(54) 



where we've dropped higher order terms that enter at 0{q'^). Thus, at leading order the influence force takes the form 



h{r) 



Sz'^ (t) 



{{SfF),), + l{{iSFF?)^), 



(55) 



The explicit form for ( 55 1 can be obtained after an integration by parts and evaluating the expectation values of the 



atom's current density separating into three distinct contributions fk = fik + f2k + fsk- 
The first component arises from the intrinsic fluctuations of the field 



hk{r) = |- dX S^'gretir, X) dk{x) {x, z"{\)) 



(56) 



zir) 



~ a 
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The retarded Green's function for the atom gret quantifies its response to an external field and so represents the 
atom's dynamic polarizability, a. The Hadamard function, , describes electric field fluctuations and is constructed 
by contracting each index of D'^ with k^. Therefore fik can be related to the recognizable form for the interaction 

energy of a polarizable body with the electromagnetic field Uint — \oiEP . 

The second component follows from dipole moment fluctuations and is the analog of the London force for the 
atom-surface geometry. Heuristically, we can construct this component of the force by considering the interaction 
energy of a quantum dipole in front of a dielectric surface. First, the atom is brought near the surface leading to 
the induction of an image dipole and therefore an image field. The interaction energy of this system is then given by 
Uint = —P' Eiraage wherc p is the atom's instantaneous dipole moment and Ei,nage is the electric field from its image. 
The image field can be expressed as the convolution of the classical electromagnetic Green's function for the electric 
field, Ql^t, with the image dipole Ef^^^g^ — J d^x' Otet p^^""^^ finally taking the expectation value of the interaction 

energy we find Uint — — J ct^x (^pj p™"^"^^ Glet- Because the atom and its image's internal degrees of freedom are 
correlated we find that the expectation value of the dipole moments of the atom and image relate directly to the 
fiuctuations of the atom's dipole moment alone 2 (^pj p™"^"^^ ~ gn — {{Qj^Qk}) showing the relation between our 
heuristic derivation and the exact expression below 



/2/c(t) = 



dX 5'=gH{r,\)dk{x)gif{x,z''{\)) 



(57) 



x—z{t) 



The last component of the force, /sfc, is similar to the first with the exception that the force arises from the induced 
fluctuations of the field. We note from the semi-classical equation of motion for the field that the electric field induced 
by fluctuations within the medium takes the form E^^^ = Jyd^x' Q'^^^ By simply taking the expression for fik 
and substituting the Hadamard function by the two-point function for induced fluctuations we find 



.2 ftf 



hkir) ^ ^ dX Sk,gret{r, X)dk{x) ({i?L(x), i5f„,(z"(A))} 



x—z{r) 



which when expanded out takes the form 



/3fc(r) = ^ r dX f d^x f dt /''di'5,et(T,A)Gr'(i,i')^;it(^"(A),a;)9,(:r')5i;'^*(2:';t',^) 
^ Jti Jv Jti Jti 

where the collective effect of the medium's fluctuations is accounted for in 



(58) 



x' —z{t) 



GT%t,t')^gH{t,t') + 



dti 



dt2 gret{t,ti)gret{t' ,t2)&H{ti,t2) 



(59) 



So we find that at leading order in perturbation theory the origin of each force component is due to the quantum and 
thermal fluctuations in a different degree of freedom. 

Equations (561, (57) and (581 which form the centerpiece of this paper, provide the fully dynamical and self- 
consistent force between an atom and a general medium including arbitrary geometry and composition |48j for the 
dielectric medium. 



VI. NONEQUILIBRIUM ATOM-SURFACE FORCE 



In this section we will study the force between an atom situated in vacuum (z > 0) and a dielectric half-space 
occupying the region [z <Q). Our goal here is to find the kernels of the electromagnetic field. The retarded Green's 
function for the field in a lossy dielectric half-space is well known, and so our focus in this section will be to calculate 
the Hadamard function quantifying the intrinsic fluctuations of the field. To do this we must find the solution for 
the time dependent field operator by solving the equation ( 44 1 . We begin by taking the divergence of the Heisenberg 

equation of motion for the homogeneous solution, Ao{s,x), for the field (44 1 giving 
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e{i,t',x)Ao{i',x) = 



(60) 



where e{t,t') is the dynamical permittivity and an over dot denotes a time derivative. 

Because the temporal gauge is incomplete the vector potential can be augmented by the gradient of any time- 
independent scalar without altering the electric or magnetic fields, the condition provided by (60) fixes the residual 
gauge freedom giving equivalent physics to the generalized Coulomb gauge [M] , 



A° = V • y dt' e{i, x)A{t', x) = 



(61) 



For a piecewise uniform permittivity the field becomes transverse within either half-space. 

For the study of nonequilibrium dynamics of our system we are restricted to compact intervals for the time in- 
tegrations appearing in the actions describing the total system. Using the Laplace transform we can simplify the 
integrodifferential equation (44) for compact time intervals which results in the transformed equation 



V^Al{s, x) -f e{is, x)s^Al{s, x) = (l)l{s, x) ^ A' {t = tj, x) + e{is, x)sA' [t ^ t^, x) 



(62) 



where e{LiJ,x) is the frequency dependent permittivity e{is,x) = 1 + Vl^pgret{is,x) and (/)o(s,x) encapsulates the effect 
of the initial conditions. We assume that the field and the medium only begin to interact at the initial time ti and 
that just prior to t = ti the field can be found in a plane wave expansion 



2(27r)3fc 



(63) 



where k = \k\, e^(fc) is the jth component of a u-polarized wave in the Coulomb gauge and aa-{k) is the destruction 
operator for a photon of wave vector k and polarization a. In the Coulomb gauge a is summed over TE (transverse 
electric) and T M ( transverse magnetic) polarization. By taking the appropriate derivatives and reading off the Fourier 
amplitudes of (63 1 we can identify (f>l,{s,k) as 



2k 



(e(is, x)s — ik)el.{k)aa-(k). 



(64) 



We now consider plane wave solutions originating from either half-space. For points within the dielectric we find 



A^-{s,k) 



e(zs)s^ -|- k^ 



(65) 



and for points in vacuum 



A^+{s,k) 



s2 -I- fc2 



(66) 



where + (-) in this section denotes waves in vacuum (dielectric) and in (66) the permittivity is taken to 1 in (?!)J(s, k). 



Inverting the Laplace transform gives the time dependence of the field amplitude for planes waves with momentum 
k. In vacuum we have 



A^+{t,k) = 



y^ei(fc)a.(fc)e-''=(*-*-)0(t~t,) 



(67) 
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and within the dielectric medium we have the general result 

A^- {t, k) = ^ f ds e^(-*-) (68) 

which depends upon the specific choice of matter-reservoir coupling. For Ohmic reservoir spectral density the disper- 
sion relation has four distinct poles and the Bromwich integral can be inverted exactly 



where ri is the residue of {e{is)s — ik)/ {e{is)s'^ + k"^) from the Zth pole and si is the /th pole. 

Given the plane wave solutions above we can now construct the total field by considering right and left-incident 
waves. Waves that approach the vacuum-dielectric interface from the right reflect and transmit. The sum of the 
incident, reflected and transmitted components gives the total solution where the boundary conditions on the field at 
the interface 

£;||(0+) = £;||(0-) and D^{Q+) ^ D^{Q-) (70) 
i/||(0+) =S||(0-) and 5,(0+) =B,(0-) (71) 

determine the Fresnel transmission and reflection coefficients. For example, in the case of TE waves the electric field 
is parallel to the vacuum-dielectric interface and because the electric field is parallel to the ve cto r potential the x and 



y components of the vector potential are continuous across z = leading to the relation (see 72 ) 1 + R^^ = T^^ (to 
avoid confusion we stress here that this is a constraint on the field amplitude and not on the energy). Because our 
dielectric model does not have a magnetic response all components of the magnetic field are continuous across z = 0. 
For the case of TE waves the continuity of the x and y components of the magnetic field give kz{l — R^e) = K-^T^^ 
where Kz is chosen so that the transmitted component of the field satisfies the wave equation in the medium. Following 
the same procedure for TM polarized waves defines the remaining reflection and transmission coefficients. This allows 
us to express the right-incident waves as 



AUx) = y Sk , 



1 



e{~k,)e^'^^'-''^+'^\y^\\a^{k)e>J{t - t 



h.c. 



(72) 



where the Fresnel coefficients become 



P-R _ k,~K, rpR _ 2k, nR _ e{k)k,-K, rpR _ 1-^ e{k)k, 

^TE — k,+K, ^TE — k,+K, ^TM — e{k)k,+K, ^TM — e{k)k,+K,- 

Here a^{k) is the destruction operator for left-moving waves with wave vector k (fc, < 0) and polarization a, the 



z-component of the wave vector in the medium is — sign(fc,)y e(fc)fc2 — k^ and the polarization vectors can be con 
veniently expressed in differential form as bte — {~^\\)~^^^i~^dy, id^, 0) and Btm — (AA||)^^/^(— —dyd^, A 



Likewise, we can follow the same procedure for left-incident waves. For a general medium one would employ (68 1, 
in the case of Ohmic reservoir spectral density we find 



ii(a;) 



a=TE,TM I 



1 



^2(27r)3fc 



0(if,)ne^'(*-*')+'^i|-^iia^(P)ei0(t-t,) 
6'(-z) + T^e'''='-^6l(z) +h.c. 



(73) 



where the Fresnel coefficients are 
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^TE — Kz+ki,z ^TE — ki,z+K, ^TM — Kz+s{isi)ki_, ^TM — e(isi)ki,z+Kz 

the wave vector in the medium is k' — (fc|| , K^) with magnitude k' = |fc'|, and fc; ^ is chosen to solve the wave equation 

in vacuum fc;.2 = ^ —si^ — fc^ with Im[fc;.2] > 0. By considering the cases where the dielectric becomes an ideal 

conductor (e -> oo), e becomes a real constant, or the dielectric becomes transparent (e = 1) one can verify that the 
expressions for the field above reduce to the expected half-space [34 ^-nd vacuum forms (after choosing an appropriate 
initial state). 

We can now construct the two-point functions of the field under the influence of the medium. The retarded Green's 



function solves (42) and can be split into two pieces; one corresponding with the source point in vacuum and the 
other with source point in the medium. This can be easily understood by noting that A^'^ is proportional to a^'^ 
and therefore (^A^A^^ = 0. Therefore we can decompose into a portion due to right and left-incident waves 

b^,,{x, x') Dli,^^{x, x') + Dl,^^{x, x'). (74) 
The retarded Green's function gives the electric field from a point dipole located at r 

Ef^x)^ I dt'g:f{x;t',r)pi{t'). (75) 



Dipole oscillations will lead to emitted radiation that will be seen at x at retarded times. By noting that the 
source point for the retarded Green's function appearing in the component of the atom-surface force due to medium 
fluctuations lies in V and the observation point lies to the future in the vacuum region we can ascertain that information 
about medium fluctuations is contained entirely in the left-incident component of ^^^^ which allows the replacement 



Qlf e[f'^ in (|58| 



The Hadamard 



mnction for the field can be constructed from the solution for the field operators ( 72 1 and ( 73 ) 



which also decomposes into a left and right-incident component 



D%{x,x') = ({A^^{x),Ai{x')}) + ({A^^{x),A^^{x')}) . (76) 

The crucial point we would like to make now is that for a dissipative medium the intrinsic fiuctuations of the field 
are exponentially suppressed in time for fluctuations emanating from within the medium. 

The time dependence of the field for wave vector k emanating from within the medium can be calculated by inverting 
the Laplace transform in (68). Because all of the poles (and/or branch cuts) in the dispersion relation lie to the left 
of the imaginary axis in the complex s-plane the intrinsic fiuctuations of the field associated with left-incident waves 
are damped away in time so that at late times 

Df',{x,x') {{Afix),Af{x')}) EE Df^''''ix,x'). (77) 
A. Long-Time Limit 

We now specialize to the case of a stationary atom and take the long-time limit. Because this system exhibits 
dissipation we expect that a steady-state exists at late times. In the previous section we have already noted that 
for general reservoir-matter coupling that field fiuctuations emanating from inside the dielectric medium will be 
exponentially suppressed in time. Within the medium the matter fluctuations will be defined by the reservoir at late 
times as the polarization field is thermalized. The collective fiuctuations of the medium is described by the kernel 
G^^'^'^{t,t') containing components due to the intrinsic fluctuations of the matter and reservoir. We showed explicitly 



for Ohmic reservoir spectral density that the intrinsic fluctuations of the matter damp away in time ( 32 ) and will for 
general reservoir-matter coupling. Therefore as — — oo we find 

/>oo poo 

GS^'^(i,t')^ / dt, dh9ret{t,ti)gret{t',h)&H{tl,t2). (78) 
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where the upper hmit of the time integrals can be taken to infinity because the matter's retarded Green's functions 
contain theta functions which restrict the integration range over past times. By taking the Fourier transform we find 
an equivalent description in frequency space 



(79) 



The retarded Green's function for the matter satisfies ( 23 1 sourced by a delta function 



^gretit, t') + U^-gretit, t' ) - j^' dt"&^'\t, t"y9ret{t" , t') = 5{t - t') 



(80) 



whose Fourier transform gives (again tf, — >■ oo) 



w'^gret(w') + U}^gret{oj') - <3'^''* {uj')gret{^') = 1 



(81) 



where gret(w') is a complex function of frequency. Multiplication of (81) by g*^^{uj') and then taking the imaginary 
part provides an identity between &^ and g^'^* [IHl [SU] 



5;et(c^')Im[(5™*(w')]5>-et(c^') = Im[5„t(^')]- 



(82) 



Using the fiuctuation-dissipation relation ( 19 ) we find 



(83) 



where the last equality holds only for ti — !■ — oo. The imaginary part of gret{^') relates directly to the imaginary part 
of the permittivity so that in the long-time limit we find the asymptotic behavior 



_,2 /-oo 



fk 



47r 



duj' S^n9:eti^')dk{x)g!'^^^{uj', f, z) + g*H{Lo')dk{x)gifitj', x, z)],=,(,) 



„2 roc 



2n 



dio' / d^x Im[£(c.')]coth(/3;fw72)g;,,(w')C:L('^'.^:a-)afe(x')C"(^',a?';x) 



(84) 



c'=z(t) 



1. Fluctuation-Dissipation Relation for the Field in Vacuum 



In the long time limit, when a steady-state has been attained, we can study the features of right-incident waves, 
, in terms of Fourier modes in the vacuum region as in (72) 



AUx) 



E 

-TE.TM 



a^(fc)e^„e^*"(*-*')+'^-^ ( 1 + R^e-''^^'' ) + h.c 



e{z). 



(85) 



The Wightman function for the field provides all information about propagation and fluctuations of right-incident 
waves in the vacuum region (evaluated here at coincident spatial arguments for later convenience) 



&^+{x-t' ,x) = (A^^{x)A'j,{t' ,x)) = 



2(27r)3 



E 



a=TE.TM 

coth l3Ek/2 COS k{t -t') sin k{t ~ t') 



(86) 



where we have assumed that the right-incident component of the field is in a thermal state at inverse temperature 
Pe- Two times the real part gives the Hadamard function and two times the imaginary part relates to the retarded 
propagator 
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&l^{x;t',x) = 2Re[&^+{x;t',x)] bl[^,^{x-,t',x) = 2ie{t-t')lm[D'^+[x-t\x)]. (87) 

After taking the Fourier transform of these stationary kernels with respect to t — t' wc find that they are related by 
the fluctuation-dissipation relation 

D'hA^'^ f , x) = 2 coth(/3£;w72) lm[bi[^ ,^{u:' , x, x)] (88) 

at the same point. 



B. Steady-State Atom-Surface Force 



For comparison with previous works we note that Q]?^^{ui\x,x') and q^g^etiy') in our treatment can be identified 
with the dyadic Green's function for the electric field, G^^ /{Air), and A'Ka{uj') the frequency-dependent polarizability 
used in [5711301 131] ■ 

By using the fluctuation-dissipation relation the atom-surface force can be brought into the same form as derived 
by others. After taking the long time limit, Fourier transforming in time, assuming global thermodynamic equilibrium 
and specifying the atom's trajectory to be static (84 1 can be used to express the force in equilibrium. This expression 
can be simplified by the use of an identity. 

If we specify that the spacetime point x lies in vacuum and to the causal future of x' with spatial component in V 



then Gl^^'^{x,x') satisfies (42) and can be shown to satisfy the identity 



d^x' ,S'^-Im[£(c.',x')]^;i;i(w',f,a;')^i:f ^(u;',f,x') = 4,ImKi,_ Ju;', f , x)] 



(89) 



For z > the fluctuation-dissipation relation requires G^'^ {oj' , z, z) = 2coth{Puj' /2)lm[Qlj*'^'{uj' , z, z)] and 



in Fourier space. 
For z > the J 

gH{<^') — 2 coth(/3a;'/2)Im[(7ret(w')]. With these relations (and borrowing the notation of others) we can rewrite fk as 

1 



fk 



duj' 5'^ coth{/3u' /2)lm[a{uj')dkix)G^j{uj' ,x, z)] 



(90) 



?(r) 



reducing it to the form for the atom-surface force as derived from MQED [50] . 

It is interesting to consider now the specific case where the temperature of the medium differs from the field. Using 
the fact that 



^^[Gret^i^', z, x)] + lTa[gl'^l'{uj\z, x)] = lTa{Ql{t{uj\ z, x)] (91) 
we can write the force in the form 



fk 



<1_ 
An 



duj' 5'^ 



2 coihiPE^ /2)g:,M)dk{x)lHQ%t{^\ z)] + g*H{J)dk{x)gif {u:' , z) 



'2tx 



dJ / d^'x Im[£(c.')](coth(/3MC^72) - co\h{fiE^' l2))gl,,{J)Q'^:, ^{J , z,x)dk{x')g:f'''{iu',x',x) 



(92) 

x — z{t) 



z(r) 



The first term is equivalent to the force associated with the atomic energy level shifts derived in [23], when the 
atom is in a thermal state at temperature Tg it gives back the atom-surface force in thermal equilibrium. The last 
term gives the correction to the force when the medium and field are out of thermodynamic equilibrium. Wc denote 
this correction to the force /I'^'' 



Jk 



2tt 



duj' / d^y 5'^gl^^{u:')dk{x) Im[e(w')] ( coth[/3M^72] - coth[/3Bw72]) 
Jv 



(93) 
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C. Nonequilibrium Steady-State Force 



In this section we elucidate the physics of the nonequihbrium correction to the force by evaluating the convolution 
of the Green's functions in ( 93 1 and comparing the results with the force in equilibrium. After some manipulation 



( 93 1 can be brought into the form 



fk {TmiTe)ew 



V2 



: ^jKe[e{uj')] - + \e{u 



duj' j dqq Lo'^Re[a{uj')] ( coth{f3M0j' /2) - coth{(3Eio' /2)) - 1 

1 , (2g2_l)(52 + |s(a;')_g2|) 



(94) 



pneqirp rp \ 



72 



X\/Re[£(cj')]-g2 + |£(a;')-g2| 



dJ dqq iu''^lm[a{uj')] ( coth(/3Mw72) - coth(^£;a;72)) ^/q^ - 1 
1 , {q^ + \e{oj')-q^\) 



lyr^ + v/el^O-g'P \e{oj')^/l^+ v/F(Z7y^|2 



(95) 



where the force splits into a component due to evanescent waves {EW) and one due to propagating waves (PW), see 
appendix for details (the properties of the polarizability have been used to simplify the frequency integrations). Note 
in particular that the evanescent component of the nonequilibrium correction is equivalent to the difference between 
two evanescent components to the atom-surface force in equilibrium (110 1. The propagating wave component gives 



rise to a constant space-independent wind term which results from the radiation of the field into the vacuum region 



(z > 0) . The imaginary part of the polarizability in ( 95 1 has support in a narrow band around the resonances of the 



atom, indeed for our oscillator model the imaginary part of the polarizability is proportional to a delta function which 
allows the frequency integral in (951 to be done directly. However, for realistic laboratory temperatures the effect of 
the propagating waves is exponentially suppressed. (See ^30] for a detailed derivation of asymptotic properties of the 
evanescent wave nonequilibrium contribution to the force.) 

We now have an expression for the atom-surface force for a stationary system when the dielectric medium is out of 
thermal equilibrium with the field. By dissecting the Lifshitz force we can intuitively argue for the form of the force 
derived previously. First, in thermal equilibrium the Lifshitz force is composed of a propagating and evanescent wave 
component. The propagating component arises from field fluctuations in the vacuum region and partially transmitted 
waves from the medium. The evanescent component arises from fluctuations of the field within the dielectric that 
partially transmit into the vacuum region from incident angles exceeding the critical angle. By subtracting the effect 
of transmitted waves from the expression for the force in equilibrium, both at temperature Te, we describe the effect 
of incident and reflected propagating waves in the vacuum region. To account for the waves partially transmitted into 
the vacuum region we need to add the effect of medium fluctuations within the dielectric at temperature Tj\/ which 
is done by adding /"'^^(Tm, Te) to the equilibrium force. This results in a modiflcation of the evanescent component 
of the equilibrium force and accounts for the blackbody radiation emitted from the surface when the global system is 
out of equilibrium 



fk = h{TE) - fr-'^'^'iTE) + f™{TM) + fr{TM,TE)pw. 



(96) 



D. Force as a function of time 



Now that we have isolated the steady-state form for the force we can identify the components giving rise to force 
dynamics. We showed previously that as time evolves the reservoir will thermalize all parties that it interacts with 
eventually resulting in a steady-state. The power of our nonequilibrium formulation is that we can follow the time- 
evolution of the force from an initial state. 

To highlight these effects, not obtainable through MQED or Lifshitz theory, we isolate the dynamical terms from 
the total force which at early times will lead to different predictions from the standard equilibrium theories. These 
corrections come in two terms. First, we showed earlier that the left- incident fluctuations of the field are damped 
away in time due to interaction with the medium. For times shorter than the inverse dissipation rate left-incident 
field fluctuations augment the equilibrium force by 



19 



Thermal Equilibrium Nor)equilibrium Steady-State Fully Nor)equilibrium 




FIG. 1. The figure above summarizes the diflterences between atom-surface interactions in equilibrium, in a nonequilibrium 
steady-state, and under fully nonequilibrium conditions. The grey bulk at left illustrates the matter where the symbols e{ij) 
and ^ denote that the medium has been coarse-grained and so its averaged effect is described by a permittivity and a stochastic 
force. The vacuum region at the right of the matter contains an atom and the field. The simplest case is thermal equilibrium 
where are all parties are described by a thermal state at temperature T, the long time limit is taken, and the position of the 
atom is fixed. The equilibrium case can be generalized by the consideration of nonequilibrium steady-states where all parties 
are described by thermal states (e.g. Ta, Tf, and Tx are the atom's, field's and reservoir's temperature), the atom is held 
fixed, and the long-time limit is taken so that the state of the entire system no longer evolves in time. A fully nonequilibrium 
treatment describes dynamics. For this case the total system is time evolved from an initial state given above as thermal for 
each party. The illustration depicts the time dependent flux of energy emanating from the medium as well as atom motion. 
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2 Jt 



(97) 
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The second contribution comes from the induced fluctuations of the field which at early times, before the matter is 
thermahzed by the reservoir, gives the additional force contribution 



fTir) = Y£'d\ dt' 5.et(r,A)5r'(^,OCt,L(^"(A),xH(a:')e:;^*'^(a:';i',x) 

An illustrated summary of the different cases treated in this paper can be found in Figure[l]. 
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Time dependent force due to field thermaUzation by eddy currents 



As an explicit example of force dynamics we'll compute the time-dependent contribution to the atom-surface force 
arising from the thermalization of the field by the medium. For this purpose we'll work explicitly with the Drude 
model which has been chosen for its simplicity. In addition we'll assume that the field is initially described by a 
thermal state, factorized from the remaining parties, where the matter is assumed to have already been thermalized 
by the reservoir. 

The Drude model dispersion relation contains two dissipative oscillating modes, the bulk plasmons, and a third 
diffusive mode arising from eddy currents [56j . For weak dissipation (dissipation rate much smaller than the plasma 
frequency) the plasmons give rise to a constant space independent force at leading order corresponding to black body 
radiation into vacuum. In distinction the eddy currents establish an evanescent field which varies with atom-surface 
spacing and so will be the focus of our investigation. After computing the Hadamard function for the field in ( 97 1 we 
find the explicit expression for the dynamical force below 



/it"(T) = -; 



^aof^W [ dX [ dk [ dqk 



coth(/3A:/2)sinf7(r- A) 
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(99) 
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where ao is the oscillator's static polarizability, we've set the initial time to zero and we've made the change of 
variables Kz = kq. In the limit of weak dissipation the root, s/, the z-component of the wave- vector, kz,l^ and the 
residue, rj, for the diffusive mode can be approximated by 



7fc2 + n] 



•Seddy ^ ~ , 2 I n2 '^eddy ~ r,n 2 i n2\ k^^eddy ^ ik^Jl - q^. (100) 



After making the change of variables y = \/\ — q^ we note that the exponential factor exp{s;r + sl,\ + i(kz,i — k^ i,)z} 
suppresses large values of y and k for Opz >> 1 at finite times. After expanding the integrand for small y and k we 
find 



/it"M«-^^£rfA [dkj\y sinn{r~X)k'cotHPk/2)y' (lOi) 

where a more precise comparison of the thermal wavelength to the time and distance scales is necessary before we can 
expand the hyperbolic cotangent. Note that the combination of exponentials in the integrand ensures convergence 

.2 1,2 

for the fc-integration and that replacing — ^2^^^ (t + A) — > + •^) results in a negligible error when >> 1. 

Making this replacement allows us to perform the y-integration which to leading order for large z gives 

/iTM^-gl^^^^^'rfA J^^dk smn{r-X)k'cotH/3k/2)e'^4^^^'\ (102) 

At room temperature (or hotter) time is the factor which determines the largest scale in the problem when r > Ins 
(for a gold surface), thereby allowing a Taylor expansion of the hyperbolic cotangent in the integrand. Subsequently 
the k and A integral can be performed where the result below is expressed to leading order for Ht >> 1 
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(103) 



A few remarks are in order. First, it is not surprising that the force depends on the inverse of the cubic distance. In 
a nonequilibrium steady-state where the field and matter are described by different temperatures it is an evanescent 
field that gives rise to the novel scaling of in the far-field [3^. The particular case we address here is the time- 
dependence of the evanescent field as it is thermalized by the eddy currents. Second, as the conductivity becomes 
very good flp/j » 1 the presence of eddy currents is suppres sed ( in this limit the Drude model approaches plasma 



model). This effect is captured by the prefactor ^^jS^p in (103|) which relates to the resistivity of the material 
indicating that the force is larger for poor conductors. 



The plots in Figure[2] quantify (1031 for a gold surface and a rubidium atom when the field is initially at room 
temperature. 



VII. CONCLUSION 



In this paper we have derived from first principles with a microscopic model the nonequilibirum force between 
an atom and a medium (modeled by a continuous lattice of non-interacting harmonic oscillators each coupled to a 
reservoir) including fully dynamical activities such as dissipation, absorption, radiation and fluctuations. In the end 
we show that the force is in agreement with the results of previous works employing MQED theory when the long-time 
limit is taken and the system has evolved to a steady state. 

Our results go beyond previous work by providing the fully dynamical force including effects of relaxation, ther- 
malization, and backaction from the field and matter in addition to being valid for general dielectric geometries and 
composition. These new results could play a role in experiments measuring the dynamical aspects of surface atom 
interactions, like quantum friction |57l 158) (and references therein) and relaxation. The dynamical effects derived 
herein are very sensitive to the model used to describe the dielectric. In particular one would expect stark differences 
between the dynamics of the force when comparing Ohmic models, like Drude, with the plasma model which neglects 
dissipation. Dynamical effects not only provide a new arena in which to explore quantum fluctuation forces but could 
aid in the resolution of controversies involving quantum friction j59H61j and the so-called thermal problem [621 163) . 
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FIG. 2. The plots above quantify the dynamical force (in Newtons) at room temperature (T = 295 K) between a rubidium 
atom and a gold surface as the field is thermalized by eddy currents. The plasma frequency and dissipation rate for gold used 
were flp = 8.9 eV and 7 = 0.0357 eV. For rubidium we set the static polarizability, ao, equal to 4.73 x 10~'^^m"^ and Q to 
rubidium's first optical resonance 2.35 x 10^^ Hz. The plot on the left shows the spatial dependence of the atom-surface force 
(neglecting the rapidly oscillating cosine term in (103l) as a function of distance in microns evaluated at different times. The 
solid line corresponds to one microsecond. From there each successive line moving upward gives the force at one microsecond 
intervals from Ifis to 4/is for the top line. The plot on the right shows the time dependence of the dynamical force as a function 
of time at an atom surface distance of one micron. The force is rapidly oscillating and so we've plotted the average force, 
indicated by the solid line, and the envelope, by the dashed lines, bracketing the force oscillations. 



As a final note we point out that beyond giving a fully dynamical description for atom motion in the presence 
of a dielectric material the nonequilibrium quantum field theoretic formulation of this problem is particularly adept 
at treating fluctuation phenomena. It will be the aim of a future study to understand the rich interplay between 
the propagating and evanescent wave components of the force and their effect upon the fluctuations of the atom's 
trajectory in space. 



VIII. APPENDIX 



A. Evanescent and Propagating Waves 



To understand the nature of the nonequilibrium force correction (93) we study the Green's function Q^l:^{uj' ,x,y) 



in the vacuum region. In equilibrium the fluctuations of the field in this region (z > 0) come in two categories; 
left-incident waves and their reflected components, and waves that partiaUy transmit into the vacuum region. By 
use of the fluctuation-dissipation relation we will use Q^I:^(io' ,x,y) to identify the effects associated with partiaUy 
transmitted flcld fluctuations. 

In the vacuum half-space Ql^*{uj\ x, y) is obtained by solving the classical equations of motion for the field, subject 
to dielectric boundary conditions in the z = plane, sourced by a delta function. Its form can be taken from [52j or 
the Appendix of [3T] and is written below for both spatial arguments lying in the vacuum half-space: 



Cy-*(L.',f,f')-e° + 



E 

a=TE,TAI 



l/c2(z+2:')gife||-(K|| 



(104) 



Above Rg- are t he Fr esnel reflection coefficients where the index a refers to the polarization (implicitly appearing 
below ( 108 ) and ( 109 )). The polarization vectors for the field in the vacuum region are e„^i{zL) where eT£;(±) = k»x z 



and eTA/(±) = ik\\z^kzk\\)/uj' , and z arc unit vectors, z normal to the vacuum dielectric interface with fc|| directed 

kn is the z-component of 



parallel to the projection of the wave vector in the plane of the interface, and k^ = uj' 

the wavevector. The free field component of the Green's functions, or bulk term, is C/° and would be present whether 
the dielectric were there or not. The bulk term leads to a divergence in the energy but has no spatial dependence at 
coincidence, and so we discard it as we require the force to vanish at infinite atom-surface spacing. 

Let us now focus on the portion of the atom-surface force in thermal equilibrium that arises from fiuctuations of 
the field, be them induced or intrinsic given by 
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(105) 



?(r) 



where FF stands for field fluctuations. Tracing over the indices of the Green's function, taking the spatial derivative 
and setting the two spatial arguments to the position of the atom we can express the FF-contribution to the force as 
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dJ / d^k\\ cj' a*(a;')coth(/3w72)Im 
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^2x n 



i2k. 



(106) 



When both spatial arguments lie in the vacuum half-space fJJ'^*(w', x, y) decomposes into two contributions. First, 
field fluctuations in the vacuum region give rise to propagating waves, in particular, waves moving in the negative 
2;-direction will reflect from the dielectric surface giving rise to waves propagating in the positive z-direction. Second, a 
field fluctuation within the dielectric producing waves propagating toward the vacuum-dielectric interface will partially 
transmit into the vacuum region. One consequence of this is that the z-component of the wave vector in ( 106 ) is not 



necessarily real. For values of |fc||| < u' we see that kz is real but becomes pure imaginary for the integration range 
|fc||| G (a;',oo). The former is associated with the propagating solutions to the wave equation in the vacuum and the 
latter with evanescent waves. 

We can isolate the influence of the evanescent modes on the force if we restrict the integration range so that the 
magnitude of the transverse momenta are strictly greater than the wave frequency 
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(107) 



where EW stands for evanescent waves. 

In this range the z-component of the wave vector is pure imaginary, kz — > in, and this property can be used to 
simplify the equation for the force by explicitly taking the imaginary part of the integrand inside the curly brackets 
{..}. We note that 
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(109) 



where k'^ — ye{uj')uj'^ — fc| is the z-component of the wave vector in the dielectric and in the last step we have used 
the identity w'^Re[e(a;')fc/] = Re[/c^](fc| + |fc^P). To bring the force due to evanescent waves into its final form we 
make the change of variables fc|| ~ uj'q and note that \/2R.e[/c^] = w Re[e(aj')]a;'^ — fc^ -t- |e(ti;')a;'^ — 
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This expression now quantifies the force on the atom due only to evanescent field fluctuations 



(110) 
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B. Nonequilibrium Correction: Relation to the Equilibrium Force 



In this section we explicitly evaluate the the nonequilibrium correction to the force beginning with the convolution 
of the Green's functions in f^^"^- Starting from (93 1 we note that the Green's functions appearing are those with 
support at one point in the vacuum half-space and the other located within the medium. Using the appendices of the 
papers [51] and [SS] we can evaluate the convolution of the Green's function product 



/C30 /•OO />U 

dx dy dzg';f''^iio',n,x)d,gZ,^{J,r2,x) (111) 
-OO J —OO J — OO 

Above Tte and Ttm are the Fresnel transmission coefficients for the TE and TM-polarized waves. After making the 
change of variables fcy = w'g, plugging in the explicit forms for the Fresnel coefficients, and setting ri and r2 to be the 
position of the atom / separates into two contributions, one from evanescent waves and the other propagating waves 
(labeled EW and PW) 
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For evanescent waves we have constrained the imaginary wave vector to be positive by expressing it as = 

Combining this with the expression for the nonequilibrium 



1 after the change of variables fcii 



U! q. 



correction to the force ( 93 1 we can express the force as ( 94 ) and ( 95 ) . 
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